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Naive dimensional analysis based on chiral effective theory, when adapted to nuclear energy 
density functionals, prescribes natural units and a hierarchy of contributions that could be used 
to constrain fits of generalized functionals. By applying these units, a large sample of Skyrme 
parametrizations is examined for naturalness, which is signaled by dimensionless coupling constants 
of order one. The bulk of the parameters are found to be natural, with an underlying scale consistent 
with other determinations. Significant deviations from unity are associated with deficiencies in the 
corresponding terms of particular functionals or with an incomplete optimization procedure. 
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Introduction. New experimental data for atomic nuclei 
throughout the nuclear mass chart are becoming avail- 
able thanks to radioactive nuclear beam efforts vi^orld- 
wide. This is coupled with increasingly sophisticated 
theoretical descriptions of low-energy nuclear phenom- 
ena [1 , '^ . These developments drive higher requirements 
for the quality and predictive power of nuclear structure 
investigations. 

Nuclear density functional theory (DFT) [3i] is the only 
available theoretical tool for the microscopic description 
of nuclear properties that spans the full nuclear mass 
chart. The non-relativistic Skyrme energy density func- 
tional (EDF) is based on local nuclear densities and cur- 
rents and is specified by a set of coupling constants. 
There are many sets of Skyrme parameters determined 
through different optimizations to experimental energies, 
radii, and other nuclear observables. After more than 
twenty years of experience, the standard Skyrme EDF 
has proved to be fairly successful in the overall descrip- 
tion of experimental data. At the same time, its limi- 
tations have become well established, and the quest for 
better accuracy and stable predictive power motivates 
going beyond the standard Skyrme functional [J, Q . 

New developments in nuclear DFT are inevitably as- 
sociated with an optimization of EDF parameters to a 
selected set of experimental data. This is problematic if 
more general density dependencies and higher powers of 
gradients lead to an explosion of new parameters without 
control over their relative importance. A possible solu- 
tion is to organize generalizations of the Skyrme func- 
tional by effective field theory principles that exploit the 
separation of scales and so establish a hierarchy of con- 
tributions. One such approach for low-energy quantum 
chromodynamics is naive dimensional analysis (NDA) for 
chiral effective field theory [g , which has been adapted to 
relativistic nuclear EDFs with encouraging results jTHlQf . 

In the NDA approach, a scaling to "natural units" is 
applied to the functional, which if successful results in di- 
mensionless parameters of order unity. This has practical 
benefits, because in numerical optimization it is advanta- 



geous to have all the parameters close to unity. But the 
use of natural units also provides guidance on maintain- 
ing a hierarchy and preventing fine tuning where higher 
orders play off against lower orders, particularly when 
linear combinations of the parameters are underdeter- 
mincd. 

The relevance of naturalness for Skyrme functionals 
was suggested long ago in Ref . [11| , but the resulting nat- 
ural units have not been widely employed (or validated) 
by DFT practitioners. The first test was very limited in 
the r ang e of functionals and focused on isoscalar parame- 
ters [ni . The goal of this paper is to investigate whether 
natural units apply more generally to existing Skyrme pa- 
rameterizations, including the modern ones, and thereby 
motivate their use in future nuclear DFT developments. 
As part of this study, an online converter to natural units 
has been created at http : //massexplorer . org, where 
one can browse and convert the Skyrme forces consid- 
ered in this paper as well as try new sets of parameters 
to check whether they are natural. 

Skyrme Energy Density Functional. The standard 
Skyrme energy density can be written in isospin represen- 
tation as a sum of kinetic and potential isoscalar (t = 0) 
and isovector (t = 1) energy density terms 

^M = |-T-o + Ho(r)-t-Hi(r), (1) 

zrn 

where the time-even part is 

Utir) = (Cfo + Cfo pD pI + C^'pt^pt + Clptn 

+ iC/ J2 + c7'ptV ■ Jt , (2) 

The isospin index t — {0, 1} labels isoscalar and isovector 
densities pt, Tt, and Jt, respectively. The standard defi- 
nitions of these local densities can be found in Ref. [1^ . 
The energy density of Eq. ^ depends on 13 parameters; 
that is, 12 coupling constants and one exponent 7, 






(3) 



which are typically obtained by adjusting the functional 
to produce certain properties in finite nuclei and/or in 
infinite nuclear matter. Historically, the Skyrme force 
and the Skyrme functional derived from it were defined 
by using the {i„, Xn} parametrization. The link between 
these two representations can be found in Ref. [3|. The 
most general Skyrme functional also contains a time-odd 
part with associated time-odd coupling constants, which 
become relevant in nuclear states with nonzero angular 
momentum (e.g., odd-mass nuclei) Q. In this work we 
consider only the time-even part and leave the time-odd 
coupling constants as a subject of future study. 

Natural units. Following Ref. [ll| , we scale the Skyrme 
coupling constants by analogy to an effective Lagrangian 
where each term is schematically written as 
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with g being a dimensionless coupling constant and /tt ~ 
93 MeV is the pion decay constant. The momentum 
A characterizes the breakdown scale of the chiral effec- 
tive theory. As such, it is expected to be in the range 
500 MeV < A < 1000 MeV. Naturalness implies that g 
should be of order unity, which in practice roughly means 
between 1/3 and 3 (unless there is a symmetry reason 
making g small). If natural, Eq. Q implies a hierarchy 
of terms with a density expansion (powers of I) and a 
gradient expansion (powers of n) |7Hll| . 

The conversion of the Skyrme couplings to natural 
units is accomplished in the present work by multiply- 
ing each by a scaling factor 



S — /•2(i-l)^n+i- 



(5) 



where I is the power of densities in the corresponding 
term and n is the number of derivatives for that term. 
In Ref. [11| only functionals with integer powers of the 
density-dependent term were considered. Here we gener- 
alize the scaling to include also fractional powers 7 used 
in Skyrme functionals by setting I — 2+j for the density- 
dependent term. At present this is just a prescription. In 
studies of dilute fermion systems in a harmonic trap, it 
was shown that terms with fractional powers in a pertur- 
bative functional followed scaling rules [ij], but this has 
not yet been derived for the nuclear case. 

We also generalize our analysis to the isovector cou- 
pling constants, which highlights the issue of possible 
additional numerical factors in the NDA prescription of 
Eqs. (|H) and ([S]). A direct extension to the isovector 
channel scales the isovector coupling constants with the 
same scale factor S as the corresponding isoscalar cou- 
plings. However, in past applications of the NDA to rel- 
ativistic meson and point coupling models, the isovec- 
tor prescription included an additional factor of four. 
This arises from the construction of the Noether cur- 
rent for an isospin transformation, which in the conven- 
tional normalization has a 1/2 with each r matrix (so 
that [Ta,Tb] = itabcTc implies Ta = Ta/2). Then the 



isovector current is l/2tp^^^Ttp and so l/2{pp — p„) is 
the corresponding charge density used in the naturalness 
analysis. Although this may be no more than a theoret- 
ical prejudice, the empirical observation in other NDA 
tests was that the scaled constants consistently came out 
closer to unity 14]. In the present study we consider 
both isovector scalings. 

More generally, to decide on possible additional numer- 
ical factors we rely on the correspondence of Skyrme EDF 
terms to those from a non-relativistic reduction of a rel- 
ativistic formulation (e.g., meson exchange with masses 
of order A) . For example, one might wonder if the spin- 
matrix u should lead to extra scaling factors between 
scalar terms and those involving the vector densities Jt- 
We find that such terms arise with the same relative fac- 
tor as terms with pr and so we scale them the same. 

For all coupling constants entering the standard func- 
tional of Eqs. ([1]) and ([2]), one has I = 2 except for the 
density-dependent constant C^j^, for which I = 2 -|- 7. 
Similarly, the power is n = for Cf, while for all 
other constants n = 2. In this way, scaling all cou- 
pling constants Cf with the associated factors S"^ , a = 
{p, Ap, r, V J} yields dimensionless constants S^CI . The 
small ranges for I and n preclude testing the fine de- 
tails of the NDA scaling hypothesis. However, by mak- 
ing a global analysis of Skyrme parametrizations, we can 
check for consistency, for trends and exceptions to natu- 
ralness, and for a preferred range of A. When parameter 
sets for extended functionals that include higher-order 
derivatives [4] and/or higher powers of density [l5| are 
available, more definitive tests of natural scaling will be 
possible. 

The list of functionals considered is given in Table ID 
In this table we also categorize the functionals based on 
the strategy used to determine the couplings. 

The test of whether we truly have natural units is 
whether S" makes the values of all scaled constants S'^C'^ 
of order unity. Their numerical values will obviously de- 
pend on the value of the cut-off parameter A [ll|. In 
our global study, the naturalness criterion can itself be 
used to extract the value of A by minimizing the devia- 
tion of the coupling constants from unity. We consider a 
logarithmic root-mean-square deviation (RMSD) 



RMSD= /l^logio|Cf(z)|, 



(6) 



because naturalness implies couplings should not be too 
small as well as not too large. If a particular coupling con- 
stant is zero, it is excluded from the logarithmic RMSD. 
In Fig. [T]we plot RMSD for 48 EDFs as a function of 
A with (scaled) and without (unsealed) the extra factor 
of four for isovector terms. It can be seen that the two 
different scalings produce different optimal A with the 
scaled result yielding a clearer minimum that is numer- 
ically more consistent with studies of relativistic func- 
tionals. However, the minima in the RMSD curves are 
quite shallow, so A cannot be considered to be sharply 
determined for the present Skyrme functionals. 



TABLE I: List of Skyrme functionals and categories used in 
this study. The categories are: a) masses of double-magic nu- 
clei (includes '"'Zr, "''Sn, ^^^Sn, and ^""Ce) used in the fit; b) 
masses of non-double-magic nuclei used in the fit; c) charge 
radii used in the fit; d) single-particle energies used in the fit; 
e) symmetric infinite nuclear matter constrains considered in 
the fit; f) asymmetric infinite nuclear matter constrains con- 
sidered in the fit; g) surface properties (neutron skin, fission 
barriers, etc.) considered in the fit; h) pairing was present 
in the fit; i) some parameters were fixed in the fit; j) pa- 
rameters extrapolated or fine-tuned from an existing force or 
functional. 
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Index 


Functionals 


Categories 


Ref 


1-2 


SkT3, SkT6 


ad i j 


[16] 


3 


SkM 


a c e f g i 


[17] 


4 


SkM* 


gj 


[18] 


5-6 


SGI, SGII 


d e j 


[19] 


7 


HFB9 


abfhi 


[20] 


8-9 


SI, SII 


a c d e f i 


[21J 


10 


SkA 


ac de g i j 


[22] 


11 


HFB16 


a b c f h i 


[23] 


12 


SkT 


a b d e g h i 


M 


13-16 


SLy4-7 


a c d e f i 


[25] 


17-18 


SkIl-2 


a b c d f g i 


[26] 


19-20 


SkI3-4 


a b c d f g i 


[26] 


21 


SkI5 


a b c d f g i 


[26] 


22-27 


MSkl-6 


abfhi 


[27] 


28-29 


SIII, SIV 


a c i 


[28] 


30-31 


SV, SVI 


a c i j 


[28] 


32-33 


SLy230a,b 


a c d e f i 


[29] 


34-39 


E, E,,, Z, Z,,, R,,, ( 


3ct a c d g i 


[30] 


40 


SkP 


a b c e f h i 


[31J 


41-42 


SkO,SkO' 


a b c d f g i 


[32] 


43 


SV-min 


a b c d g h 


[33] 


44 


SkOy" 


i j 


M 


45 


SkMP 


aj 


[35] 


46-47 


SkX, SkXc 


a b c d e f 


[36] 


48 


RATP 


a d e f i 


[37] 



In the present study, we choose to use the scaled 
isovector coupling constants for which the optimum is 
A = 687 MeV (but the precise value does not affect our 
conclusions). In Fig. [2] we have plotted the scaled cou- 
pling constants for all the functionals of Table HI Also, 
we plot the square-roots of individual RMSD contribu- 
tions given by the functional to the total RMSD value. 
It can be seen that the Skyrme functionals have almost 
all of their parameter values in the interval (1/3, 3) with 
the bulk between 1/2 and 2. Exceptions are discussed 
below. 

We also make a comparison between different represen- 
tations of two particular functionals: SIII and IIFB16. 
The parameters of these functionals are listed in Ta- 
ble |IT] first by using the (i, a;)-parametrization and then 
by the natural units parametrization, obtained from the 
corresponding coupling constants. As can be seen, in 
the (t, a;)-parametrization these two functionals seem to 
be quite different from each other. However, when ex- 
pressed in natural units the coupling constants of SIII 
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FIG. 1: Logarithmic RMSD as a function of A with (scaled) 
and without (unsealed) an extra factor of four for isovector 
terms. See text for details. 



and HFB16 are order unity. In Table HIl we also list in 
natural units the average, minimum, and maximum value 
for each coupling constant found in the set of 48 function- 
als. This information may provide useful insights into the 
expected values and ranges of coupling constants for fu- 
ture attempts to fit new functionals. 

Deviations from order unity. The deviations of the 
coupling constants C^ from unity are illustrated in the 
summary plot in Fig. [5] As noted earlier, almost all pa- 
rameters are found to lie within the interval of (1/3,3). 
In terms of naturalness, we do not observe any signifi- 
cant differences between the functionals that are strictly 
based on the Skyrme force and the extended function- 
als. However, significant deviations still exist for some 
particular Skyrme functionals for the coupling constants 
Cq,Ci '' , and C^'\ and more generally for C^-^, which 
appears borderline unnaturally large in many cases. If we 
accept that nuclear functionals are characterized by nat- 
uralness, such deviations could indicate some deficiency 
in the associated term of the functional or they could sim- 
ply reflect a specific strategy applied to the optimization 
procedure. 

While in some cases no definite cause has been identi- 
fied, we can identify various examples of unnatural cou- 
plings that do have probable explanations. For example, 
it is well known that tensor terms are rather poorly con- 
strained by the experimental data; most Skyrme func- 
tionals do not include tensor terms at all. The signif- 
icant deviations seen in Fig. [2] for the Cq parameters 
most likely refiects this situation. 

Another instructive example is the deviation for C-^ '' 
in the case of Skll (case A in Fig. [5]) . The optimization 
of Skll excluded the isotope shift data while all other 
Ski forces (SkI2-SkI5) consider these data. This isotopic 
shift data contains the charge radii difference between 
4°Ca and ^^^Ca, and ^"Spb and ^i^Pb. All Ski functfonals 
are, however, optimized by using diffraction radii data, 
which is closely related to charge radii. 
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FIG. 2: (color online) Scaled coupling constants \Ct\ at A =687 MeV (top) and contributions of individual functionals to the 
total RMS value (bottom). The filled symbols refer to the isoscalar coupling constants and empty symbols to the isovector 
ones. The ordering of functionals by index is the same as in Table U 

TABLE II: Comparison between the SIII and HFB16 functionals in the (f , a;)-parametrization and using natural units. We 
also list for each coupling constant the average, minimum, and maximum value found in the set of 48 functionals. In all cases, 
A = 687 MeV was used. 





{t, x)-parameters 






Couplings 


m natural units 








SIII 


HFB16 




SIII 


HFB16 


Average 


Min. 


Max. 


to 


-1128.75 


-1837.23 




-0.4767 


-0.7759 


-0.7977 


-1.2380 


-0.4465 


tl 


395.0 


383.521 


fip 


1.2076 


1.9295 


1.7656 


-0.9795 


4.5761 


t2 


-95.0 


-3.41736 


fip 


0.7623 


0.7509 


0.7824 


0.2723 


1.2616 


ti 


14000.0 


11523.0 


'-'ID 


-3.0493 


-2.3825 


-2.2010 


-6.0116 


1.9790 


Xo 


0.45 


0.432600 


C(T 


0.6059 


0.4464 


0.5606 


-0.0856 


2.9421 


Xl 





-0.824106 


CI 


-1.6726 


-0.2048 


-0.3626 


-2.9469 


3.5160 


X2 





44.6520 


f-iAp 


-0.8597 


-0.8702 


-0.8765 


-1.7491 


-0.4636 


X3 





0.689797 


cf" 


0.9301 


-0.8998 


-0.5531 


15.5202 


2.5762 


Wo 


120.0 


141.100 


(^0 


-1.2288 


-1.4449 


-1.2385 


-1.6384 


-0.8957 








c7' 


-1.6384 


-1.9265 


-1.0875 


-2.1846 


5.4272 








ci 





1.1300 


0.5159 


-0.6021 


1.4212 








ci 





1.3208 


1.6502 


-5.0026 


3.6049 



The deviation in Cg for SV results from an artificially 
imposed vanishing density-dependent term, which results 
in a too-low isoscalar effective mass (0.38). The EDF 
RATP demonstrates a clear example of an anomalously 
small Cf^ « -0.0019 (not seen on the scale of Fig. [2|). 
This can probably be attributed to the fitting procedure, 
where the focus was mainly on infinite nuclear matter 
properties for astrophysical applications. Similarly, C^ ^ 
of SkMP (case C) is also very small. This can be linked 
to the fact that this functional was obtained by mixing t- 
and x-parameters of SkM* and SkP functional and mak- 
ing small adjustments to the volume part of the func- 
tionals. Therefore, almost no attention was paid to the 



surface part either in RATP or SkMP functionals. 

Another example is an anomalously small Cg, found 
in SkX and SkXj. (cases D and E). This is due to the 
fact that in fitting these forces, much emphasis was put 
on the single-particle energies. This leads to an effec- 
tive mass close to one, and therefore a small coupling 
constant. Similarly, in the MSk series (case B) the fit 
favored effective mass close to one, and it was therefore 
set by hand either to 1.0 or 1.05. This, however, does 
not imply that single-particle energies are not suitable 
observables in the fitting procedure. 

Finally, the seemingly unnaturally large CCp couplings 
may reflect an inadequate treatment of terms with frac- 
tional density dependence. Alternatively, the fact that 



the isovector Cfp couplings are also sometimes unnatural 
hints at a problem with the scaling of terms associated 
with pion-range physics. The density matrix expansion 
applied to the leading long-range contributions from chi- 
ral effective field theory may shed light on this issue. 

These examples illustrate that the use of natural units 
in nuclear DFT not only introduces the simplicity of di- 
mensionless coupling constants and the convenience of 
their order-unity values, but also can give valuable point- 
ers to potential deficiencies of the physics invoked when 
constructing and optimizing the functional. 

Conclusions. In this study, the coupling constants of 
a large set of Skyrme EDFs have been examined for nat- 
uralness as an extension of Ref. [ll[. While the lim- 
ited range of density and gradient terms in the standard 
Skyrme parameterizations means that a definitive test of 
NDA chiral naturalness is not possible, the best function- 
al are consistent with naturalness and a scale A of about 
700 MeV. Significant deviations from unity can be asso- 
ciated with deficiencies in fitting the functional or with 
specific optimization procedures. This motivates using 



naturalness as a guiding principle for constructing new 
generalized Skyrme functionals. An online natural units 



converter has been set up at http : //massexplorer . org 
as a tool for such applications. Further investigation is 
needed to better understand how to treat non-analytic 
density dependence and hybrid functionals where the 
density matrix expansion is used for long-range contri- 
butions. Finally, we note that the phenomenologically 
successful finite-range Gogny functionals can be accom- 
modated within the Skyrme framework [38|, which can 
be used to broaden the application of natural units. 
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